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SUMMARY 

The FORTRAN 7094 computer program for the calculation of the design-point perfor- 
mance of turbojet and turbofan engine cycles is presented along with the thermodynamic 
equations used. This program requires as input the airplane Mach number, the altitude- 
state conditions, turbine-inlet temperature, afterburner temperature, duct -burner tem- 
perature, bypass ratio, coolant flow, component efficiencies, and component pressure 
ratios. The output yields specific thrust, specific fuel consumption, engine efficiency, 
and several component temperatures and pressures. The thermodynamic properties of 
the gas are expressed as functions of temperature and fuel to air ratio. Sample input and 
output are included for an example case. 


INTRODUCTION 

Advanced aircraft for supersonic flight, high-payload long-range subsonic flight, and 
vertical flight require the development of advanced propulsion systems. The study of 
system characteristics and components for such airbreathing propulsion systems has 
been undertaken at the Lewis Research Center. As part of these studies, thermodynamic 
analyses of design-point performance of turbojet and turbofan engine cycles must be per- 
formed. To assist in these analyses, a computer program was developed. 

The program requires as input the airplane Mach number, the altitude -state condi- 
tions, turbine-inlet temperature, afterburner temperature, duct -burner temperature, 
bypass ratio, coolant flow, component efficiencies, and component pressure ratios. The 
thermodynamic properties used in this analysis are expressed as functions of tempera- 
ture and fuel to air ratio. The fuel is assumed to be of the form (CHg) . The results of 

the analysis include specific thrust, specific fuel consumption, engine efficiency, and 
several component temperatures and pressures. 


The equations used in the analysis and the FORTRAN 7094 computer program are , 
presented in this report along with a description of the input and output. Sample input and 
output are included for an example case. 


METHOD OF ANALYSIS 

An analysis for the performance of a general engine is derived herein. The perfor- 
mance of this engine can be studied in terms of specific fuel consumption and specific 
thrust. The general engine consists of components of turbojet and turbofan engines. An 
analysis of the performance of these specific engines is obtained from the general engine 
analysis by proper specification of the input. 


Cycle Description 

The general engine cycle is shown in figure 1. The air enters the diffuser and the 
major portion of its velocity head is changed into a pressure head. The lower velocity 
air then enters the fan and is compressed. The air flow is then divided into the main 
flow and the duct flow. The main flow enters the compressor and is further compressed. 
The major portion of the main flow then enters the combustor with fuel, and combustion 
occurs. The small remaining portion of the main flow bypasses the combustor and is 
used to cool the turbine. The combustor exit gases are then expanded in the turbine, 
producing work to drive the compressor and fan. The turbine exit gases and the coolant 
flow then mix and enter the afterburner with fuel, and combustion occurs. The hot gases 
are then expanded in the main nozzle to produce thrust. The duct flow enters the duct- 
burner with fuel, and combustion occurs. The hot duct gases are then expanded in a noz- 
zle to produce additional thrust. 



Figure 1. - Schematic diagram of general engine. 
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Derivation of Equations 


The equations used in the analysis of the general engine are derived in this section. 
The thermodynamic properties used for this analysis are functions of temperature and 
fuel to air ratio. The specific heat of the gas is expressed by a polynomial equation. 
Appropriate integrations of this equation yield the enthalpy change and the entropy func- 
tion. The entropy function is herein defined as 


A <p = 



(All symbols are defined in appendix A. ) The derivation of the equations for specific 
heat, enthalpy change, and entropy function are given in appendix B. Since the specific 
heat is a function of temperature and fuel to air ratio, it is expressed as C p (T, f) in this 
analysis. Since 


and 



C p (T,f)dT 



Cp(T, f) 

T 


dT 


they are expressed as Al^Tg, Tj,f) and AcpCTg, Tj,f), respectively. When the fuel to 
air ratio is zero, these quantities will appear with temperatures only. The fuel is as- 
sumed to be of the form • 

Engine inlet . - The static inlet conditions of the diffuser T Q and p Q are a function 
of the altitude. The inlet velocity is 


V 0 * M 0 Vg VT 0 


( 1 ) 


where Mq is the Mach number at which the airplane is traveling. The specific heat 
ratio is 


3 


y 


l 



1 " C p< T 0> J 


The total temperature at the inlet is 


T, o= T 0 + 



where 


C 


P 


Ah < T 0> t q) 

T i- T o 


'( 2 ) 


(3a) 


(3b) 


The correct total temperature is then obtained by an iterative procedure involving equa- 
tions (1) to (3b). The total pressure at the diffuser inlet is obtained from 


R. 


— 1" rr = AvdJ, T 0 ) 


P T 


Therefore, 


Av(r 0 ,T 0 )J/R a 

0 ~ 0 

Diffuser . - Since there is adiabatic flow in the diffuser, 


( 4 ) 


Ti = Tq (5) 

The pressure ratio across the diffuser r j q is an input parameter, which is a function 
of the Mach number. Therefore, 



( 6 ) 


An example variation of this parameter is presented in reference 1. 

Fan . - The fan pressure ratio P'^t/P’^ is a variable. To determine the ideal fan- 
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exit temperature, the isentropic flow equation is used. Therefore, 


R„ 


Pi- 


A^ F = ~ In — 7- 
b J Pj 


and 


A<p(Tj*^ T^) - A<pp 

Therefore, Tjt id can be determined from equation (7b). The fan work is 

_Ah( T » 1 , >id , T i) 


Ah p = 


(7a) 


(7b) 


(8a) 


where fan efficiency is a design parameter. The fan-exit temperature T'^i is deter 
mined from 


Ah(Ti*, T'j) = Ah p 


(8b) 


The total airflow is 


W tot = W 


a, D + 


W 


a, m 


(9a) 


where W r, is the duct airflow and W„ „ is the main airflow. The ratio of the duct 
a, D a, m 

airflow to the main airflow W„ r ./W 0 „ is called the bypass ratio b. Therefore, 

a, ±j a, m 


w 


tot 


(1 + b)W 


a, m 


(9b) 


Compressor. - The overall pressure ratio of the fan and compressor is a variable. 
Thus, the compressor pressure ratio is 



V 


jyn 

pi./pi 


( 10 ) 


The ideal compressor-exit temperature Tl, ic j is obtained from 


A «’ (T 2, id 




( 11 ) 
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The compressor work is 


Ah c = 


_ Ah ^ T 2, id’ T l'^ 


(12a) 


where the compressor efficiency is a design parameter. The compressor exit tempera- 
ture Tl, is determined from 


Ah(T^, Tj») = Ah c 

Combustor. - An energy balance for the combustor yields 


(12b) 


W f,m h f 


+ 77 B W f HVF 


+ (W 


a, m 


W c> h a 


(W 


a, m 


w„ + w f )h 


f/“g 


(13) 


For the enthalpy of the fuel to be equal to zero, the enthalpy reference temperature T^ 
must be equal to the temperature of the incoming fuel. As discussed in appendix B, the 
enthalpy change of the gas can be expressed as 


Ah = (Ah + Ah.f) -L- 
g a D 1 + f 


(14) 


Therefore, substituting equation (14) into equation (13) and dividing by (W - W ) 

^ lit c 

yield 


W f HVF W, 

w ’ w + ih < T 2' T R> - Ah < T 3> V + w ih b< T 3' T R> 

w a, m w c w a, m w c 


The fuel to air ratio based on combustor airflow is 


W 


f, m 


W - W 
a, m c 


Ah(T^, T r ) - Ah(T^, T r ) 
HVF, b - Ah^, T r ) 


(15) 


where the turbine inlet temperature Tg and the combustor efficiency are design 
parameters. The fuel to air ratio based on main flow is 


W 


f , m 


W 


f, m 


/W - W 
' a, m c’ 


W 


f, m 


W„ W - W 
a, m a, m c 


W„ 


a, m 


W„ 


- W, 



(16) 
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where the coolant -flow ratio W /W o is a design parameter. The combustor pres- 

t a,y m 

sure ratio r^ ^ * s an input parameter. 

Turbine. - The turbine work is equal to the fan work plus the compressor work: 


(W 


a, m 


- W c + W f, m> Ah T - W tot Ah F + W a, m Ah C 


(17) 


Substituting equation (16) into equation (17) and solving for turbine enthalpy drop yield 


Since 


Ahrp = 


(1 + b)Ahp + Ah^ 


W 


1 - 


W 


a, m 



W 


f, m 


W - W 
a, m c ; 


Ah T' 3 ,T’ 4 , 


W 


f, m 


W - W 
a, m c/ 


1= Ahn 


(18a) 


(18b) 


the turbine exit temperature T^ is determined from equation (18b). To determine the 
turbine pressure ratio, the ideal turbine-exit temperature T 4. id is needed. Therefore, 
the ideal turbine work is 


Ah 


T , id 


Ahrp 

v T 


and 


T 3’ T 4,id’ 


ly m 

W - W 
w a,m w cy 


= Ah, 


T, id 


(19a) 


(19b) 


where turbine efficiency is a design parameter. The ideal turbine-exit temperature 
T 4 id is determined from equation (19b). Thus, the turbine pressure ratio is 


— = exp 

P i 


A^( T 3, T 4,id’ 


W. 


f, m 


\M„ — I 


W - W 
a, m c/ 


g 


(R 


(20) 


where the equation for the molecular weight of the gas is given in appendix B. 
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Coolant mixing . - The turbine-exit gas and coolant-flow mix just downstream of the 
turbine. This mixing is assumed to take place without any change in the mainstream 
total pressure. However, there is a change in total temperature. An energy balance for 
the mixing section with 0° R as the reference temperature yields 


W c Ah(T',0) + (W am -W c + W f> m )Ah fa 0, 


w < m 
t, m 

W - W 
a, m 


/ W 

< W a,m +W f,m) ah Km.°. ^ 
\ a, my 


Dividing by W & m yields 

/ W„ 
c - Ah^, 0) + 1 1 - c 


W 


a, m 


W 


a, my 


(i + ' VJt o, Wf ' m 

lw -w/\ w - w / 

\ a, m c/ \ a, m c/ 


/ W, \ / W, 

a, m/ \ a, my 


Solving for the exit enthalpy yields 


“’-•■a-- 


w„ / \v 

c Ah(T2,0)+ll 


W 


a, my 


1 + W «.° 

l w -w / \ w -w 

\ a, m c/ \ a, m c 


W, m 
1 + f ’ — 
W 

a, m 


( 21 ) 


Total exit temperature T^ M is determined from equation (21). 
The total pressure is 



( 22 ) 


Afterburner . - Three cases are considered for the afterburner. For case I with no 
afterburner, 
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(23) 


P 5 = P 4 T 5 = T 4M 


W 




w 


a, m 


For case n with the afterburner not lighted, 


5,4,n a 5 “ a 4M 


W 


f,AB_ Q 


W 


(24) 


a, m 


where r^ ^ n is given. For case IH with afterburning, an energy balance on the after- 
burner yields 


W. 


W. 


W, 


W r 


^AB 


HVF + ( 1 + ]h’,,„ = (l + _Ji£L + 1 ’ AB 


W 


a, m 


W„ 


4M 


a, my 


W Q m W 0 m 
a, m a, m, 


(25) 


Solving equation (25) for the afterburner fuel to air ratio yields 


W / W \ / W \ 

Ah(T’ T R ) + Ah.(T’ , Tj-lii” - 1 + -AiH ]Ah T’ „, T R , -L™. ) 

t, AB . 5 R 5 Rw a, m \ w a.m/ V 4M R W a ,J 


w 

w 


a, m 


»AB HVF - 4h b< T 5' T R> 


(26) 


where the afterburner temperature Tg is a design parameter. The total-pressure ratio 
across the afterburner is 


P’ c 


5 _ r 

r - r 5,4 


(27) 


which is a design parameter. Therefore, the total fuel to air ratio of the mainstream is 


W f,tot_ W f,A B| W f,m 
W a,m W a,m W a,m 


(28) 


Main nozzle. - Full expansion is assumed in the mainstream nozzle. Therefore, 
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( 29 ) 



and 


Since 


A<? n = -JL- In 


JM 


g 


Pc 


A <p 




W f,tot\ 
W a, m / 


= a<p n 


(30a) 


(30b) 


the ideal nozzle exit temperature Tg can be determined from equation (30b). The 
main nozzle -exit velocity is 



(31) 


where i// is the velocity coefficient and a function of the airplane Mach number (ref. 1). 

The equations derived thus far are for the main flow. The equations for the duct flow 
are now derived. 

Duct burner. - Three cases are considered for the duct burner. For case I with no 
duct burner, 


P? = P'^ Tlj = T*i» 


W 


f,D _ 


= 0 


W 


a, D 


(32) 


For case n with the duct burner not lighted, 


p^7 = r 7» l'» n T 7 " T l r 


W 


W 


f,D - 

a, D 


(33) 


For case in with duct burning, an energy balance on the duct burner is 

W a, D ah < T i’ * T R> + W I, DB<’'DB>< HVF > * ("a, D + W f , DB> Ah 8 


(34) 
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Dividing equation (34) by W„ n and solving for the duct -burner fuel to air ratio yield 

• cl. XJ 


W f , DB _ Ah ( T 7> t r) - Ah ( T i t > t r) 

r?DB HVF - A V T 7> T R) 


(35) 


where the duct -burner temperature is a design parameter. The duct -burner fuel to 
air ratio based on the main flow is 


W f, DB _ W f , DB W a, D _ W f , DB b 


W„ 


a, m 


w a,D w a,m w a,D 


The total pressure ratio is 


P 7 

P i' 


- r 7, 1’ 


(36) 


(37) 


which is a design parameter. 

Duct nozzle. - Full expansion is also assumed in the duct nozzle. Therefore, 


and 


Since 



(38) 


<R 


a< ^d n — =■ 

u ’ w JM, 



(39a) 


^8, id’ T 7’ N 


(39b) 


the ideal exit temperature Tg ^ is determined from equation (39b). The duct nozzle- 
exit velocity is 


v 8 ~ ^ 1/ 2gJ Ah ( T 7’ T 8, id’ 


W 


f, DB 


W 


(40) 


a, D 
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where ip is a function of Mq (ref. 1). 

Specific thrust. - The specific thrust of the engine is defined as the net thrust divided 
by the total airflow: 


ST = 


( W 8 V 8 + W 6 V 6 - W tot V 0> 
gW tot 


(41) 


Substituting the weight -flow relations yields 


ST = 


("a. D + w l, DB> V 8 + <»a. m + W f . tot> v 6 ~ < x ♦ b > w a, m v 0 


g(l + b)W, 


a, m 


(42) 


and rearranging equation (42) yields 


ST = 


[ b + Mv 8+ 1 + ^]v 6 ., 1+ b,v r 

a, m / \ a, m 

g(l + b) 


(43) 


Specific fuel consumption . - Specific fuel consumption is defined as the total fuel 
flow in pounds per hour divided by the net thrust in pounds: 


SFC = 


Wf.tot + Wf.DB) 360 ^ 


("a, D + W f , DB> V 8 + <»», m ♦ W f , tot> V 6 ' <» + b > W a, m V 0 


(44) 


Dividing equation (44) by W Q yields 

Cly 111 


SFC = 


+ w t.PB\ 3600g 
V W a,m W a, m. 


' b + V 8 + f 1 + + b,V ' 

a,m/ \ a, m, 


(45) 


Engine efficiency . - Another performance parameter that is often used is engine effi- 
ciency: 
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(46a) 


Thrust power 
Heat added 


The thrust power is equal to the net thrust multiplied by the inlet velocity. The heat 
added is 


W f (HVF)J 

Therefore, 

3600V 0 

Ve (SFC)(HVF)J 

The method presented for determining the performance of a general jet engine ap- 
plies to several engines: dry turbojet, afterburning turbojet, dry tuibofan, and duct- 
burning turbofan. The performance of any one of these engines is obtained by eliminating 
the appropriate components of the general engine from the calculation procedure. 

The analysis of a dry turbojet is obtained by eliminating the duct equations, setting 
the fan pressure ratio (Pj»/Pj) equal to 1, the bypass ratio (b) equal to 0, the fan effi- 
ciency (t? f ) equal to 1. 0, and taking case I for the afterburner section. The analyses of 
the afterburning turbojet are the same as the dry turbojet except that case Hor HI for the 
afterburner section is used. The analyses of the turbofan engines are obtained by elimi- 
nating the afterburner, that is, by setting P' 5 = P^ M , and Tg = T| M . 

DESCRIPTION OF INPUT 

A description of the input for the FORTRAN 7094 computer program is given in this 
section. The input quantities are turbine-inlet temperature, coolant -flow ratio, overall 
pressure ratio, fan pressure ratio, turbine efficiency, compressor efficiency, combustor 
efficiency, afterburner efficiency, duct-burner efficiency, bypass ratio, afterburner tem- 
perature, duct -burner temperature, airplane Mach number, altitude, temperature at this 
altitude, pressure at this altitude, combustor pressure ratio, afterburner pressure ratio, 
duct -burner pressure ratio, velocity coefficient of nozzle, diffuser pressure ratio, and 
heating value of fuel. The input is arranged so that parametric variation of the design 
parameters can be run. 

The input format with sample data is shown in table I. 


(46b) 


(47) 
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TABLE I. - Concluded. INPUT FORMAT FOR FORTRAN 7094 COMPUTER PROGRAM 
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The definitions of the terms used on the cards follow: 


Card 1: 

R32 

R54 

R54N 

ETAC1 

ETAT1 

HVF 

NT3 

NP2 

Card 2: 

WCWA1 

EAB1 

EB1 

TR 

NM 

Card 3: 

EDI 

DED 

EDM 

R711 

R711N 

NP11 

NB 

Card 4: 

DETAC 

ETACM 

DETAT 

ETATM 

DWCWA 

WCWAM 

DEB 

EBM 

DEAB 

EABM 


combustor pressure ratio 
afterburner pressure ratio when afterburning 
afterburner pressure ratio when afterburner is not lighted 
compressor efficiency, first value or only value 
turbine efficiency, first value or only value 
heating value of fuel, Btu/lb 

number of turbine-inlet temperatures (must be an integer) 
number of overall pressure ratios (must be an integer) 


coolant -flow ratio, first value or only value 
afterburner efficiency, first value or only value 
combustor efficiency, first value or only value 
enthalpy reference temperature, °R 

number of airplane Mach numbers, must be equal to number of altitudes 


duct -burner efficiency, first value or only value 
increment of duct -burner efficiency 
maximum value of duct -burner efficiency 
duct -burner pressure ratio when burning 
duct -burner pressure ratio when not burning 
number of fan pressure ratios, must be an integer 
number of bypass ratios, must be an integer 


increment of compressor efficiency 
maximum compressor efficiency 
increment of turbine efficiency 
maximum turbine efficiency 
increment of coolant -flow ratio 
maximum coolant -flow ratio 
increment of combustor efficiency 
maximum combustor efficiency 
increment of afterburner efficiency 
maximum afterburner efficiency 
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Card 5: 

4 

ETAF1 

DETAF 

ETAFM 

Card 6: 
AM 

Card 7: 
VCI 

Card 8: 
R10A 

Card 9: 
ALTI 

Card 10: 
TOSA 

Card 11: 
POSA 

Card 12: 
P2P1 

Card 13: 
P11P1 

Card 14: 
BA 

Card 15: 
T3 

Card 16: 

T71 

DT7 

T7M 

T51 


fan efficiency, first value or only value 
increment of fan efficiency 
maximum fan efficiency 


airplane Mach number, maximum of 50 


nozzle -velocity coefficient, one for each Mach number 


diffuser pressure ratio, one for each Mach number 


altitude, ft, maximum of 50 


temperature of air at given altitude, °R 


pressure of air at given altitude, inches of mercury 


overall pressure ratio, maximum of 50 


fan-pressure ratio, maximum of 50 


bypass ratio, maximum of 50 

turbine -inlet temperature, °R (maximum of 50) 

duct -burner temperature, °R (first value or only value) 
increment of duct -burner temperature, °R 
maximum duct -burner temperature, °R 
afterburner temperature, °R (first value or only value) 


DT5 increment of afterburner temperature, °R 

T5M maximum afterburner temperature, °R 

KODE determines engine type, must be an integer 

LD index on input statements, must be an integer (If LD = 1, all 16 data cards 

are read for the next set. If LD = 2, the next set of data consists only 
of card 16. ) 

KT5 determines when afterburner is used, must be an integer 

KT7 determines when duct burner is used, must be an integer 

The values of KODE, KT5, T51, KT7, and T71 for the various engines are given in 
table n. 


TABLE II. - VALUES OF KODE, KT5, T51, KT7, AND T71 


Engine 

KODE 

KT5 

T51 

KT7 

T71 

Dry turbojet 

0 

0 

0.0 

0 

0.0 

Afterburning turbojet, afterburner 
not lighted 

0 

-1 

0.0 

0 

0.0 

Afterburning turbojet 

0 

2 

t ab 

0 

0.0 

Dry turbofan 

2 

0 

0.0 

0 

0.0 

Duct -burning turbofan 

2 

0 

0.0 

2 

t db 

Duct -burning turbofan, duct burner 
not lighted 

2 

0 

0.0 

-1 

0.0 
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DESCRIPTION OF OUTPUT 


An example of the output obtained is given in table HI. The first four lines are input 
data, except for Tl, total temperature at fan inlet and PI, total pressure at fan inlet. 

The next line contains the titles for the following columns: 

P2P1 overall pressure ratio, which is an input and primary parametric variable 

WFWA overall fuel to air ratio, (lb fuel)/(lb air) 

P4P3 turbine total-pressure ratio 

ST specific thrust, (lb thrust)/(lb air) 

SFC specific fuel consumption, (lb fuel)/(lb thrust)(hr) 

Til fan-exit or compressor -inlet total temperature, °R 

T2 compressor -exit total temperature, °R 

T4 turbine-exit total temperature, °R 

T5 afterburner temperature (If there is no afterburning, T5 is equal to T4M, the 

total temperature of the gases at the exit of the mixing station. ), °R 
T7 duct -burner temperature (If there is no duct burning, T7 = Til. ), °R 

Pll fan-exit total pressure, psia 

P2 compressor-exit total pressure, psia 

P3 combustor-exit total pressure, psia 

P4 turbine-exit total pressure, psia 

P7 duct nozzle-inlet total pressure, psia 

V6 main nozzle-exit velocity, ft/sec 

V8 duct nozzle-exit velocity, ft/sec 

ETA-E overall engine efficiency 
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TABLE m. - SAMPLE OUTPUT 
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PROGRAM VARIABLES 


ALT 

altitude 

AMI 

airplane Mach number 

B 

bypass ratio 

CP(T, F) 

specific heat, function of temperature and fuel to air ratio 

CPI 

temporary storage of CP(T, F) 

DDH 

difference between Ah and Ah(T 2 , T^,f) 

DDS 

difference between A <p and A<p(T 2 , T^,f) 

DDSC 

difference between A q> and A<p(T 2 , T^, f) for compressor 

DHC 

enthalpy change across compressor 

DHF 

enthalpy change across fan 

DHT 

enthalpy change across turbine 

DH4M 

enthalpy at station 4M 

DSC 

eq. (11) 

DSF 

eq. (7) 

DSN 

eq. (30) 

EAB 

afterburner efficiency 

EB 

combustor efficiency 

ED 

duct-burner efficiency 

ETAC 

compressor efficiency 

ETAF 

fan efficiency 

ETAT 

turbine efficiency 

GAMMA 

ratio of specific heats 

h(t 2 , t v F) 

enthalpy change of gas 

ha(t 2 , Tj) 

enthalpy change of air 

HF(T 2 , T t ) 

Al^ of eqs. (14) and (B8c) 

POPOS 

total- to static-pressure ratio at diffuser inlet 

POS 

pressure of air at flight altitude 

P2P11 

compressor pressure ratio 


P5P4 

pressure ratio across afterburner 

P7P11 

pressure ratio across duct burner 

RIO 

temporary storage of R10A 

8(T 2 ,T v F) 

entropy function 

TOS 

temperature of air at flight altitude 

TO 

total temperature at diffuser inlet 

TOl 

temporary storage of TO 

Till 

another name for Til 

T11D 

temporary storage of Til 

T11ID 

ideal fan-exit total temperature 

T21 

temporary storage of T2 

T21D 

ideal total temperature of compressor exit 

T21D1 

temporary storage of T21D 

T31 

temporary storage of T3 

T41 

temporary storage of T4 

T41D 

ideal turbine-exit total temperature 

T41D1 

temporary storage of T41D 

T4M 

total temperature at station 4M 

T4M1 

temporary storage of T4M 

T5 

nozzle -inlet total temperature 

T6S 

static temperature at main nozzle exit 

T6S1 

temporary storage of T6S 

T7 

duct nozzle-inlet total temperature 

T8S 

static temperature at duct nozzle exit 

T8S1 

temporary storage of T8S 

VO 

airplane velocity 

VC 

nozzle -velocity coefficient 

W(F) 

molecular weight, function of fuel to air ratio 

WCWA 

coolant flow ratio 
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WFAWA 

WFDWA 

WFMWA 

WFWA1 

WFWA 


fuel to air ratio of afterburner 
fuel to air ratio of duct burner 

fuel to air ratio of main combustor, based on compressor airflow 
fuel to air ratio in combustor 
total fuel to air ratio 
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PROGRAM LISTING 


JET ENGINE CYCLE PROGRAM 

DRY TURBOJET KODE = 0 T5= 0.0 KT5 = 0 
TURBOJET AFTERBURNER NOT LIT KODE=0 T5=0.0 KT5=-l 
AFTERBURNING TURBOJET KOOE = 0 T5 = AFTERBURNER TEMP KT5 = 2 
DRY TURBOFAN KODE = 2 T7 = 0.0 KT5 = 0 KT7 = 0 

DUCT-BURNING TURBOFAN KOOE = 2 T7 = D.B. TEMP KT7 = 2 KT 5 = 0 
TURBOFAN DUCT BURNER NOT LIT K0DE=2 T7=0.0 KT7=-l KT5=0 

DIMENSION T3! 53 ) * P 2P H 50 ) . T5A ( 5 0) , AM 150) , VC I (53 ) , RIO A 1 50 ) , 
ITOSAI 50),P0SA{ 50). ALT I ( 50 ) ,P1 1 PI ( 5 0) , BA ( 50) 

C THERMODYNAMIC PROPERTIES 

CPI T» F ) =! 1 . 0/( i.0-t-F))*lAC-BB*T+CC*T**2-D0*T**3 
l«-EE*T**4+( AAA+BRB*T-CCC*T**2+DDD*T**3 
2-EEE*T**4)*F) 

W ! F ) = ! i.O+FJ/I .034522+. 035648*F) 

H!T2«T1»F)=! I.O/I 1.0+F) )*! AC*! T2-T1)- IBB/2.0) *!T2**2 
l-Tl**2)+!CC/3.0)*l T2**3-T1**3)-!DD/4.0)*1 T2**4- 
2T 1**4 ) + t EE/5.0)*! T2** 5- Tl**5 ) + (AAA*!T2-T1) 

3+1 BBB/2. ) *( T2**2-Tl**2)-(CCC/3.0)*( T2**3-T1**3) 

4*1 DDD/4.0 )*( T2**4-Ti**4)-(EEE/5.0)*(T2**5-Tl**5) ) 

5*F ) 

S(T2*Tl»F)=(1.0/( l • 0+F ) ) * ( AC*A LOG ( T2 / T1 ) — BB* 

KT2-T1 )+<CC/2.0)*( T2**2-TI**2)-(DD/3.0)*(T2**3-T1**3) 

2*1 EE/4.0)*! T2**4-T1**4)+F*!AAA*ALUG!T2/T1)+ 
3BBB*!T2-Tl)-!CCC/2.0)*IT2**2-Tl**2)+(DDD/3.0>* 

4! T 2**3-T 1**3 )- 1 EEE /4 .0 )* ( T2**4-T1 * *4) ) ) 

HA! T 2* T 1 ) =AC* ( T2- T 1 )- ( RB / 2 . 0 ) * ( T2**2-T1 **2 ) 
l+ICC/3.0) *(T2**3-T1**3)-!DD/4.0)*!T2**4-T1**4) 

2+! EE/5.0)*! T2**5-T1**5) 

HF! T2« T 1 ) = ! AAA* ! T2-T1 ) +( BBB /2. 0) * ( T2**2-T1 **2 ) 

1-! CCC/3.0)*! T2**3-Tl**3)+!DDD/4.0)*! T2**4-T1**4) 

?-! EEE/5. 0 ) *( T2**5-T1**5) ) 

AC=. 24062 
BB= .0 1 7724E- 3 
CC=.038506E— 6 
DD= .0 12662E- 9 
EE=.0013012E-12 
AAA=. 22091 
BBB=. 51 822E-3 
CCC=. 19462E— 6 
DD0=.O45089E-9 
EEE* .0043 275E- 12 

10 READ! 5*50)R32»R54,R54N,ETAC1»ETAT1 ,H VF , N T3 , NP2 ,NT5 

50 FORM AT (6F7. 3*313) 

RFAD! 5. 54)WCWA1*EAB1.EB1»TR*NM 

54 FORMAT! 4F6.1. 13) 

READ! 5. 56)ED1,0ED.EDM,R711,R711N,NPU *NB 

56 FORMAT! 5F6. 1,213) 

READ! 5,55) DETAC , E TACM ,DE TAT ,£ TATM ,DWCWA ,WCWAM ,D£B , E 8M , DE AB, EABM 

55 FORMAT! 10F6.1 ) 

READ! 5, 53)ETAFl,0ETAF,ETAFM 
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READC 5*53) t AM( I )* I=l*NMJ 
. RFAD( 5* 53 ) { VC 1 ( I )* I = l*NM) 

READC 5*53)(R10A< I )* I =1*NM) 

REAO( 5* 53)C ALTIC I )*I=1,NM) 

READC 5* 53 ) ( TO SA ( I )*I = 1*NM> 

REAO( 5* 53 ) C POSA ( I )* I = 1*NM) 

53 FORMAT { 12F6.3) 

READC5*52)(P2P1(J)* J = 1*NP2) 

52 FORM ATC 12F6.2) 

READC 5* 52 MP11P1C J> *J=1*NP11) 

READC 5*51)CBA(K)*K =l*NB) 

READC 5* 51HT3C J )* J=1*NT3) 

51 FORMAT ( 12F6.1) 

9 READC 5*57)T71*0T7*T7M*T51 *D T5 * T5M * KODE *LD*KT5*KT7 
57 FORM AT ( 6F6.1* 413) 

R=53.3 
CJ -778 • 0 
G=32.2 

00 103 I = l.NM 
TOS = TOSAC I) 

POS = POSAC n 
ALT = ALT IC I) 

AMl=AMC I ) 

VC * VC1I I > 

RIO = R10AI I) 

ETAC = ETAC1 
ETAT = FT ATI 
EB = EB1 
WCW A = WCWA1 
EAB = EAB 1 
ED ^ EDI 
ETAF = ETAF1 
T5 =* T51 
T 7 = T 71 
VCD = VC 

POS = POS* .49115 
400 IF ( KODE .EQ. 0) GO TO 500 
IF CKT71 64*65*66 

65 WRITE! 6* 303 ) 

303 FORMAT C 1H1* 20X* 33HDR Y TURBOFAN ENGINE CYCLE PROGRAM) 

P7P11 = 1.0 

GO TO 80 
64 WRITEC6.304) 

304 FORMATC 1H1 * 25X* 29HTURB0FA N ENGINE CYCLE PROGRAM) 

P7P11 = R 71 IN 

GO TO 80 

66 WRITEC 6* 305) 

305 FORMATC 1H1* 20X* 41HDUC TBURNI NG TURBOFAN ENGINE CYCLE PROGRAM) 
P7P11 = R711 

GO TO 80 
500 WRI TEC 6* 200 ) 

200 FORMATC lHi* 30X * 29HTURB0 JE T ENGINE CYCLE PROGRAM) 

IF CKT5) 63*62*61 
62 WRI T EC 6* 300 ) 

300 FORMATC 1H0* 14HN0 A FTERBURNER * 5X* 1 1 HP5 /P4 = 1.0) 

GO TO 80 


61 WR I T F C 6, 301 ) R54 

301 FORMAT! 1H0. 17HW ITH A FTERB URNI NG . 5X , 7H P5 / P4 =,F5.4) 
GO TO 80 

63 WRITE! 6, 302) R 54N 

302 FORMAT! 1HG, 35HNO AFTERBURNING WITH AN AFTERBURNER, 
15X.7HP5/P4 =, F 5.4 ) 

80 CONTINUE 

IF!KT5.F0. 0) P5P4 = 1.0 
IF! KT5. EQ. ! -1 )) P5P4 = R54N 
IF IKT5.GT. 0) P5P4 = R54 
WRITE! 6, 206) 

206 FORMAT! 1H0, 3X, 1HM, 10X, 3HALT) 

WRITE! 6, 207 ) A Ml, ALT 

207 FORMAT! 1H,F6. 3, F12.1) 

IF ! AM 1 .EO. 0.0) GO TO 30 
GAMMA= 1.0/C 1 . 0-R / ! CP ! TO S , 0. 0 ) *C J ) ) 

VO = AM1*S0RT( G*R*GAMMA*TOS) 

CPl=CP! TOS, 0.0) 

TO 1 = 0.0 

2 T0=T0S+V0**2/( 2.0*G*CJ*CP1) 

CPl = H!TO,T0S,0.0)/! TO- TOS ) 

IF ! ABS ( TO-TO 1 ) .LT. .5) GO TO 5 
T01 = TO 
GO TO 2 

5 POPOS=EXP! S!T0.T0S,0.0)*CJ/R) 

GO TO 40 
30 TO =T0 S 
VO = 0.0 
POPOS = 1.0 
C DIFFUSER 
40 T 1 = TO 

PO = P0P0S*P0S 
PI = R10*PO 
WRITE! 6, 1010) 


1010 

FORMAT ( 

117HJ 

TOS 

ETAC 

ETAT 

ET AF 

HVF 

1 RIO 

R 32 P 7 1 1 

VC 

VCD 

WCWA 

VO 

Tl 

PI 


2 POS ) 

WRITE! 6. 1020) TOS.ETAC ,E TAT ,ET AF ,H VF , Ri 0 , R32 , R7 1 1 , VC ,VCO,WCWA, 
lVO.Tl.Pl.POS 

1020 FORMAT ! 1H , F 7. 1 , 3 I F 7. 3 ) ,F 9. 0 , 5 ! F 7. 3 ) , F 8 . 4 , 2 l F8 . 1 ) , F8 . 1 , 

1F7.3) 

WRITE! 6, 199) 

199 FORMAT! 1H0 , 2X , 2HEB , 5X , 3HEAB , 4X , 2HE 0 ,4 X ,4HP5 P4 ) 

WRIT F ( 6, 201 )EB» EAB * ED * P5P4 
201 FORMAT! IH , F6.3.3F7.3) 

81 DO 101 K = 1 ,NT3 
T31 = T3(K ) 

DO 102 L = 1.NP11 
C FAN 

OSF = ( R /CJ )4 ALOG ( P 1 IP 1 ( L ) ) 

T11ID1 = 0.0 

T11ID = T1*EXP!DSF/. 24062) 

310 DOS = OSF-S! T1 1 ID.T1.0.0) 

T 1 1 ID = T11ID + DnS*TllID/CPCTlllD,0.0) 

IFCABStTllID— T11ID1) . L T . . 5 1 GO TO 311 
TllIOL = T11ID 
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311 


GO TO 310 

OHF = HIT11ID.T1.0.0)/ETAF 
T 1 10 = 0.0 

Til = OHF/. 24062 + Ti 
312 OOH = OHF - H( T ll.Tl «0.0) 

Til = Til + DOH/CP( Til ,0.0) 

IF (ABS( T11-T110).LT..5) GO TO 314 
T 1 1 0 = Til 
GO TO 312 
314 Till = Til 

DO 104 KB = l.NB 
9 = BA(KR) 

WRITE) 6. 1030) 

1030 FORMAT! 36HJ B PUP 1 T3 ) 

WRITE! 6, 1040) B , P 1 IP l ( L ) , T3 ( K ) 

1040 FORMAT! 2XF6. 3, F9. 2.F9. 1) 

WRITE! 6, 1050) 

1050 FORMAT! 132H P2P1 WFWA P4P3 ST SFC Til T2 T 

14 T5 T7 Pll P2 P3 P4 P7 V6 V8 

2 ETA-E ) 

DO 100 J = l» NP 2 
C COMPRESSOR 
T21 = 0.0 
T21D1 = 0.0 

P2P 1 1 = P2P 1! J ) /P 1 IP 1 ! L ) 

DSC = !R/CJ)*AL0G!P2P11) 

T21D = Tl l*EXP! DSC/. 24062) 

18 OOSC » DSC-S!T210»Tll»0.0) 

T2 1D=T210+DDSC*T2 10/CP l T2 10,0. 0) 

IF! ABS! T21D-T21DI ) .LT. . 5) GO TO 6 

T2101=T21D 

GO TO 18 

6 DHC = HI T210.T11.0.0) / ETAC 
T 2=DHC/ . 24062+Tll 

110 DDH=DHC-H!T2 .Til, 0.0) 

T2=T2+0DH/CP( T2,0.0) 

IF! ASSi T2-T21 ) .LT. .5) GO TO 7 
T21=T2 
GO TO 110 
C COMBUSTOR 

7 WFW Al=! HA! T31.TR )-HA I T2, TR ) ) / l H VF*EB 
l-HF! T31.TR ) ) 

WFMWA=WFWA1*I l.O-WCWA ) 

C TURBINE 

OHT =( DHC +< l.+B )*DHF ) /! ( 1.+WFWA1)* (l.-WCWA) ) 

T 4=T 3! K )-DHT/. 24062 
T41 = 0.0 

15 OOH = H! T31.T4.WFWA1 ) - OHT 
T 4=T 4+ OOH /CP ( T4.WFWA1) 

IF! ABS! T4-T41) .LT. .5) GO TO 16 
T41 =T4 

GO TO 15 

16 OHT 10= DHT /ETA T 

T 41D = T3l K )-OHT ID/. 24062 
T 41 01 = 0.0 

112 OOH = H! T31.T410.WFWA1) - DHT10 
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T41D=T41D+DDH/CP( T41D.WF WA 1 ) 

IF( A BS I T 41D-T41D1 ) .LT . .5) GO TO 111 
T 41D1=T 410 
GO TO 112 

111 IF ( T41D) 1007, 1007, 1111 

1111 P4P3=EXPIS(T410.T31.WFWA1)*WIWFWA1)*.50346) 

IF ( WCW A .EO. 0.0) GO TO 24 
C COOLANT MIXING 

DH4M=( WCW A*H( T2 , 0. 0. 0. 0)t ( l.-WCWA) * ( 1 . + WF WA1 ) * 
1H1T4.0.0.WFWA1) )/< l.+WFMWA ) 

381 T4M=DH4M/. 24062 
T4MI = 0.0 

25 D0H=DH4M-H( T4M , 0. 0, WFMWA ) 

T4M=T 4M+DDH/CP ( T4M. WFMWA ) 

IFIABS(T4M-T4M1).LT. .5) GO TO 26 
T4M 1 =T 4M 

GO TO 25 
24 T4M = T 4 

26 IF ( KT5.NE . 0 ) GO TO 17 
C AFTERBURNER 

T 5 = T4M 
P5P4 = 1.0 
WFW A = WFMWA 
WAWF=l . O/WFMWA 
GO TO 19 
17 CONTINUE 

IF I KT5.NE.I-l) ) GO TO 20 
T 5 = T4M 
P5P4 = R54N 
WAWF=1. O/WFMWA 
WFW A = WFMWA 
GO TO 19 
20 CONTINUE 

IF ( T5 .LT. T4M) GO TO 99 
P5P4 = R 54 

W FAW A= ( HA ( T 5, TR ) +HF { T5 , TR ) * WF M WA- ( 1 . + WF MWA ) 

1 *H( T4M, TR. WFMWA ) ) /( HVF*EA8-HF ( T5.TR) ) 

WFW A=WFAWA+WFMWA 
WAWF= 1 . O/WFWA 
19 Pll = P 1*P11P 1(L ) 

P2 = P2P11*P11 
P 3 = P 2*R 32 
P 4 = P4P3*P3 
P5 = P5P4*P4 
C MAIN NOZZLE 

P0SP5 = P0S/P5 
P6S = POS 

IF ( P 5 .LT. P6S ) GO TO 1005 
IFiWFWA.GT. .06763) GO TO 105 
DSN= 1 . 9862*AL0G ( P6S/P5) /W( WF WA ) 

T6S=T5*EXP( OSN/.24062) 

T6S l = 0.0 

11 DDS = 0SN-S1 T6S, T5.WFWA) 

T6S=T6S+DDS*T6S/CP( T6S.WF WA ) 

IF CABS! T 6S-T6S1 ) .LT..5) GO TO 12 
T6S 1 = T6S 
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GO TO 11 

* 12 V6=VC*SQRT( 2 • 0*G* C J* H ( T5 * T6S • WF WA ) ) 

IF ( KODF.EO. 0) GO TO 414 
C DUCT-BURNER 

319 IF { KT7.NE. 0) GO TO 320 

T 7 = Till 
P7P11 = 1.0 
WFDWA = 0.0 
GO TO 350 

320 IF ( KT7.NE.(-in GO TO 321 
T 7 = Till 

P7P11 = R 71 IN 
WFDWA = 0.0 
GO TO 350 

321 IF { T7.LT.TIU GO TO 98 
P7P11 = R 71 1 

WFDWA = (HA! T7» TR )-HA( T 111* TR ) ) / ( H VF*ED-HF ( T7 » TR ) ) 

C DUCT NOZZLE 
350 P7 = P 7P 1 1*P 1 1 

DSN = 1.9862*AL0G!PGS/P7) /W1WFDWA) 

T8S = T7*EXP!DSN/. 24062) 

T8S1 = 0.0 

411 DOS = DSN-SI T8St T7* WFDWA ) 

T8S=T8S*DDS*T8S/CP< TBS* WFDWA) 

IF { ABS! T8S-T8SI) .LT..5) GO TO 412 
T8S1 = T8S 
GO TO 411 

412 V 8 = VCD* SORT ( 2 . 0*G*C J *H ( T7. T8S * WFDWA ) ) 

WFDWA - WFDWA*B 

ST = I ! l.+WFMWA)*V6+!8*WFDWA)*V8~(l.+8)*V0)/<G*(i.+B) ) 

WFW A = WFDWA ♦ WFMWA 

SFC = 3600 .0*G* WF WA/ ! f 1 • +WF MWA )*V6+(3+WFDWA) *V8-(1.*B)*V0) 

WFW A = WFWA/( 1 .*B ) 

GO TO 415 

414 ST = ( { 1.0 + l.Q/WAWF )*V6-V0) /G 

SFC = 3600. 0*G*( 1.0/WAWF > /! <1. 0+1. 0/ WA WF > *V6- VO) 

415 ETAE = 3600.0* VO /I SFC*H VF*C J ) 

IF ( KODE .NE. 0 ) GO TO 416 
V8 = 0.0 

P7 = 0.0 

T 7 = 0.0 

416 WRITE! 6 t 1060) P 2P 1 < J ) ♦ WF WA ,P4P 3 t ST ♦ SFC * Tl 1 ♦ T2 , T4 , T5 , T7*Pll.P2 f 
1P3#P4*P7*V6»V8*ETAE 

1060 FORMAT! 1 X. F5. 2 ♦ F8. 5 . F 7.4.F 7. 2 , F 8. 4 *5F 7. 1 , 5F 7. 1.2 F7 . 1 ,F6 . 3 ) 

GO TO 100 
99 WRITE! 6* 210) 

210 FORMAT! lH0tl9HT5 IS LESS THAN T4M) 

GO TO 100 

98 WR I TE(6»2ll> 

211 FORMAT! 1H0. 19HT7 IS LESS THAN Til) 

GO TO 100 

105 WRITE! 6. 106) T3IK), P2P11J),WFWA 
106 FORMAT! 1H0. 4HT3 =* F8. 1 . 3X • 7HP2 PI = ,F4. 1 ,4X , 7HWF/W A =,F8.6) 

GO TO 100 

1005 WR I T E ( 6 » 1006) 

1006 FORMAT! 1H0*21H P5 IS LESS THAN P6S ) 


GO TO 100 

1007 WRITE! 6# 1008) 

1008 FORMAT! 1H0, 18H T41D IS NEGATIVE ) 
100 CONTINUE 

104 CONTINUE 
102 CONTINUE 
101 CONTINUE 

T 5 = T 5 + 0T5 
IF ( T5.LT.T5M) GO TO 81 
T5 = T 51 
17 * T7 ♦ DT7 
IF ( T7.LT. T7M) GO TO 81 
IF (ED.GE.EDM) GO TO 405 
£0 = EO + DED 
GO TO 400 

405 IF < WCWA .GE. WCWAM I GO TO 401 
WCW A = WCWA ♦ OWC WA 

GO TO 400 

401 IF ( EB .GE. EBM ) GO TO 402 
ER = EB + DER 

GO TO 400 

402 IF [ FAB .GE. EABM ) GO TO 403 
EAB = EAB ♦ OEAB 

GO TO 400 

403 IF (ETAT .GE. ETATM) GO TO 404 
ETAT = ETAT +DETAT 

GO TO 400 

404 IF ( ETAC.GE. ETACM) GO TO 406 

ETAC = ETAC + GET AC 
GO TO 400 

406 IF! ETAF.GE.ETAFM) GO TO 103 
ETAF = ETAF+OETAF 

GO TO 400 
103 CONTINUE 

GO TO ( 1 0* 9) ♦ L D 
END 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, November 22, 1966, 
720-03-01-35-22. 
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APPENDIX A 


b 

C P 

f 

g 

HVF 

h 

Ah 

Ah b 

J 

M 

M 

P 

R 

<R 


r 

SFC 

ST 

T 

V 

W 

W 


SYMBOLS 


bypass ratio 

specific heat, Btu/(lb)(°R) 

fuel to air ratio, (lb fuel)/(lb air) 

gravitational constant, 

32.17 ft/sec 2 

heating value of fuel, Btu/lb 

enthalpy, Btu/lb 

enthalpy change, Btu/lb 

correction to air enthalpy 
(eq. (B8c)), Btu/lb 

mechanical equivalent of heat, 

778 ft-lb/Btu 

Mach number 

molecular weight 

pressure, psia 

gas constant, ft-lb/(lb)(°R) 

universal gas constant, 
ft-lb/(lb mole)(°R) 

pressure ratio 

specific fuel consumption, 

(lb fuel)/(lb thrust)(hr) 

specific thrust, (lb thrust)/(lb air) 

temperature, °R 

velocity, ft/sec 

weight flow, lb/hr 

weight of component per pound 
of air, Ib/lb air 


y ratio of specific heats 
tj efficiency 

A <p entropy function 

i p velocity coefficient 

Subscripts: 

AB afterburner 

a air 

B combustor 

C compressor 

c coolant 

D duct 

DB duct burner 

e engine 

F fan 

f fuel 

g gas 

id ideal 

m main 

N nozzle 

n not lighted 

R reference 

T turbine 

t thrust 

tot total 

0 diffuser inlet 


1 


afterburner outlet 
nozzle outlet 


fan inlet 
l' fan outlet 

2 compressor outlet, main burner 

inlet 

3 main burner outlet, turbine inlet 

4 turbine outlet 

4M mixing station 


5 

6 

7 duct -burner outlet 

8 duct nozzle outlet 
Superscript: 

(') total, as applied to state points 
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APPENDIX B 


DERIVATION OF COMBUSTION GAS THERMODYNAMIC-PROPERTY EQUATIONS 

Reaction Stoichiometry 

The fuel used was of the form C n H 2 n . Therefore, the general combustion equation 


is 


C n H 2n + - n0 2 - nC0 2 + nH 2° 
Eliminating n, the reaction and the formula weights are 


(Bl) 


CH 2 + ^° 2 - co 2 + H 2 0 
14?026 48^00 44?0K) 18?016 


(B2) 


For f pounds of fuel used, the amount of Og used is 

48.000 


14.026 


f = 3.422 f 


The amount of CC >2 formed is 


44 --P 10 f = 3. 138 f 
14.026 


and the amount of HgO formed is 


1JL.Q16 f = 1.284 f 
14.026 

The weights of the components of air per pound of air are oxygen, 0. 2314; nitrogen, 
0.7552; argon, 0.0129; and carbon dioxide, 0.0005. Therefore, the net amounts of com- 
ponents left after the reaction of f pounds of fuel with 1 pound of air are 

W n = 0.2314 - 3.422 f 
u 2 
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= 0.0005 + 3. 138 f 



W H 2 0 =L284 f 


W M = 0.7552 
W 2 


W Ar = °' 0129 


and the weight of the gas is 1 + f pounds per pound of air. 


Specific Heat 


The specific heat of the gas is 



Z( wc p), 


components 
1 + f 


(B3) 


The specific heat of each component is expressed in the form (ref. 2) 


Cp = A + B(Txl0" 


} ) + c(txio“ 3 ) 2 + d(txio -3 ) 3 + e(txio" 3 ) 4 


(B4) 


where Cp is in Btu per pound mass per °R and T is in °R. The coefficients A, B, C, 
D, and E, are obtained from reference 2. These coefficients and the molecular weight 
of each component are given in table IV. 

TABLE IV. - COEFFICIENTS AND MOLECULAR WEIGHTS OF COMPONENTS 


Components 

Coefficients in eq. (B4) 

Molecular 

weight, 

M 

A 

B 

C 

D 

E 

Nitrogen 

0.25410 

-0.032690 

0.047685 

-0.014973 

0.0014885 

28.016 

Oxygen 

. 20334 

. 029680 

.0089971 

-.0058842 

. 00076764 

32.000 

Carbon dioxide 

. 11097 

.21110 

-.088140 

.018003 

-.0014317 

44.010 

Water vapor 

. 44266 

-.033155 

.087761 

-.024552 

.0021734 

18.016 

Argon 

. 12438 

0 

0 

0 

0 

39.944 
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, Therefore, the equation for the specific heat of the combustion products is obtained 
by substituting equation (B4) into equation (B3). The resulting equation is 

I 2 3 

(c p ) = Vo. 24062 -0.017724 Tx10" 3 + 0.038056x(io' 3 t) - 0. 012662x(io' 3 t) 


+ 0. 0013012x1 


:(i0' 3 t) 4 + [O. 

- 0. 19462x(lO' 3 T) 2 + 0.045089 x(io" 3 t) J - 0 . 0043275x(i0" 3 t) 


22091 + 0. 51822X10 U T 
\3 


f> (B5) 


C p (T , f) =(C p ) 


Enthalpy Change 

The enthalpy change can be expressed as 

*T 0 


Ah 


-f 

■/T. 


CL dT 


Substituting equation (B5) into equation (B6) and integrating yield 


a V 


-_i.ro. 24062(T 2 - Tj) - °-°H724xlO- 3 . T 2) + 0, 


038056X10 


-6 


(B6) 


(i - t D 


0.012662X10' 9 ^,4 _ T 4^ + 0.0013012X10 12 ^,5 _ T 5j + 22091(T 2 - Tj) 

0.51822X1Q- 3 ^,2 _ t 2^ _ 0. 19462X10' 6 ^,3 _ + 0^045089X10^ ^,4 _ 

2 ^ 4 

0.0043275X10 *“ ( B 7) 


- 


( T t - T l) 

f ► 


* 


5 


The enthalpy change of the gas can also be expressed as 


Ah g = < Ah a + Ah b f) (B8a) 

where Ah & is the enthalpy change of 1 pound of air 

Ah a = Ah(T 2 , Tj) (B8b) 

and Ah^ is the correction to the air enthalpy due to combustion and is expressed as 


Ah fe = 0.22091(T 2 - Tj) + — 5 1§, 22X1 ° 3 (t| 

2 


T 2j _ 0. 19462X1Q" 6 ^3 
3 



° ^°89 X 10- 9 ( t 4 . T «) . 0. . 0043276X10- 12 ^5 . T S) (B8c) 
^ 5 


Entropy Function 


For isentropic flow, 



(B9) 


This equation is used in the evaluation of ideal processes in turbomachines. For con- 
venience, the right side of equation (B9) will herein be called the entropy function. The 
entropy function is 



Substituting equation (B5) into equation (BIO) and integrating yield 


(BIO) 
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A<p tT = —L- Jo. 24062 In — - 0. 017724xlO' 3 (T„ - TJ + 
g 1 + f 1 T x 1 1 


0.038056X10' 6 U2 m 2 


( T 2 - T D 


0.012662X10 -9 ^,3 _ + (L_0013012xi0" 12 / m 4 m 4> 


( T 2 - T t) 


+ f 


0. 22091 In — + 0. 51822X10" 3 (T 2 - Tj) - °- 19462X10 6 ^2 _ 


-TN 


0. 045089xl0~ 9 4p3 _ T 3\ _ 0, 0043275X10~ 12 
3 V 2 17" 4 


( T 2 - T l) 


} (Bll) 


A(^(T 2 , T 1 ,f) = A<p 


Molecular Weight 

The molecular weight of the gas is equal to the weight of the gas divided by the total 
number of moles (sum of the moles of the components). Therefore, molecular weight 
can be expressed as 

M — (B12a) 

y P) 

/ ^ /components 

The resulting equation is 


M = L±i (B12b) 

g 0.034522 + 0.035648 f 
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